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Abstract 

A new model for a spin 1/2 ladder system with two legs is introduced. It is 
demonstrated that this model is solvable via the Bethe ansatz method for arbitrary 
values of the rung coupling J. This is achieved by a suitable mapping from the 
Hubbard model with appropriate twisted boundary conditions. We determine that 
a phase transition between gapped and gapless spin excitations occurs at the critical 
value J c = 1/2 of the rung coupling. 
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Research in spin ladder systems continues to attract considerable attention. This is 
primarily motivated by the desire to understand the phenomenon of high temperature 
superconductivity observed in doped antiferromagnetic materials. In studying ladder 
materials, it is anticipated that knowledge of two-dimensional systems can be gained 
by way of extrapolation from the one-dimensional scenario, where there exists a greater 
understanding of the physics from the theoretical perspective. Moreover, it is possible to 
experimentally study ladder materials and numerical simulations are easier to treat than 
for the two-dimensional case which facilitates a greater interaction between theory and 
phenomenology. For a review of these aspects we refer to |TJ . 

More recently, there has been a push to apply the mathematically rich techniques of the 
Quantum Inverse Scattering Method (QISM) B and associated Bethe ansatz procedures, 
which have successfully been used in the study of one-dimensional quantum systems Q, 
to provide some insight into the behaviour of the ladder systems. In order to extend 
the standard one-dimensional approach of the QISM to the case of ladders, a number of 
approaches have thus far been proposed. 

Recall that central to the QISM is the Yang-Baxter equation, solutions of which, known 
as i?-matrices, provide a systematic means to produce a quantum system with an infinite 
number of conserved charges which guarantees integrability. In the algebraic formulation, 
there is an associated algebra known as a Yangian and in this approach the Bethe ansatz 
is described in terms of the representation theory of this algebra. The work of Frahm 
and Rodenbeck Q developed a construction for ladder systems where the extension from 
the one-dimensional system to the ladder was obtained by an algebra homomorphism of 
the Yangian algebra. In this manner, the symmetry algebra of the ladder system remains 
the same as the original one- dimensional model. Closely related to this approach is that 
adopted by Muramoto and Takahashi || who employed the higher order conservation 
laws obtained through the QISM to define a two-leg system which generalizes the long 
celebrated Majumdar- Ghosh zigzag ladder model ||. Such zigzag systems can also be 
obtained via an alternating spin formulation such as that considered by Abad and Rios 
0] for the t — J model. 

Alternatively, the approach can be considered where the symmetry algebra is extended 
to describe the ladder model. This notion was promoted by Wang M who successfully 
constructed a two- leg ladder system based on the symmetry algebra su(4) as opposed to 
the su(2) symmetry of the one-dimensional Heisenberg chain. A great advantage in em- 
ploying this method is that it allows for the introduction of rung interactions of arbitrary 
coupling strength which in turn allows for the identification of the phase diagram in terms 
of this parameter. The rung interaction is introduced by way of a chemical potential like 
term. In Wang's analysis, two critical values of the rung interaction were found indicating 
three phases, one of which exhibited gapped spin excitations. Subsequently, this method 
was extended by Batchelor and Maslen || who determined a class of n-leg systems each 
with a corresponding gl(2 n ) symmetry and has also been generalized to some other cases 
0,0,0. 

Yet another direction was taken by Albeverio et. al. ]TS| who searched for an i?-matrix 
solution which was not previously associated with a known one- dimensional system. They 
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were able to determine such a solution where the rung interactions appear in the Hamil- 
tonian directly from its derivation from the i?-matrix. They could then generalize this 
system to one with two free interaction parameters similarly to Wang's solution discussed 
above. A subsequent analysis of this model has been undertaken in WA . 



Our aim in this work is to demonstrate that it is possible to obtain a solvable ladder 
system with arbitrarily coupled rung interactions directly from an i?-matrix solution. To 
achieve this end, we begin with the coupled spin formulation of the Hubbard model as 
introduced by Shastry [|Tj|, on a closed lattice with twisted boundary conditions. The 
algebraic Bethe ansatz solution of this model has been studied by Martins and Ramos 
T6fl . By means of carefully chosen transformations, we map this model to a spin ladder 



system with periodic boundary conditions. Remarkably, the resulting model assumes a 
simple form with three basic forms of interaction. The energy expression in terms of a 
Bethe ansatz solution is also obtained. For this model, the rung interactions are not of 
the chemical potential type referred to above. Rather, the rung interaction parameter 
appears explicitly in the Bethe ansatz equations and so it is reasonable to expect that 
the behaviour of this model to differ from the class of ladder models with chemical po- 
tential type rung interaction. We find the critical value J c = 1/2 for the rung interaction 
parameter indicating the transition between gapped and gapless phases. 

Here we show solvability of the following two-leg ladder Hamiltonian with an even 
number of rungs and periodic boundary conditions. Explicitly, the global Hamiltonian is 
of the form 



L-l 



H = ^ + h L1 (1) 

i=l 

where the local Hamiltonians read 

= (a+o-j + *r4) ( r ^y +1 + ( T tT + T * r i) ( a ^y + j / 2 + • ( 2 ) 

Above, the coupling J can take arbitrary values. 

The ladder system is depicted graphically below. Across the rungs there is the usual 
Heisenberg XXX interaction while along the legs the interactions alternate between pure 
and correlated XX exchanges. Clearly the correlated exchange is a four body interaction 
depending on the spins of the opposing leg. In the J = limit the correlated exchanges 
have no real physical significance since for this case the model may be mapped back to 
two decoupled XX (or free fermion) chains with twisted boundary conditions. This is 
in some contrast to the case of M where in the absence of rung interactions the model 
maintains non-trival interactions between the legs. In the thermodynamic limit, the 
boundary conditions become irrelevant and we conclude that this region is gapless. On 
the other hand for large J the situation is the same as the two-leg Heisenberg ladder. In 
this limit the ground state consists of a product of rung singlets and the excitations are 
gapped Hence we expect there to exist a finite critical value of J defining the phase 
transition. 
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Ti-l 



A more detailed analysis of the model can be made using the fact that there exists an 
exact solution. The energy levels of this model take the form 

N 

E = AJN -3JL + J2 2 cos k j ( 3 ) 
where the variables kj are solutions of the following Bethe ansatz equations 

3 = 1,2, ...,N 
Z = 1,2,..,M. (4) 

The states associated with solutions of the above equations are eigenstates of the total 
spin operator 

with eigenvalues — 2M. 

The existence of the critical point is evident from the Bethe ansatz equations. Exact 
diagonalization of the two-site Hamiltonian shows that there is a unique ground state 
for J > 1/2 which is given by the product of the two rung singlets. For L sites it then 



- (-1) N exp(iLkj) 



M 

n 



sin kj — ui + iJ 



sin kj — ui — iJ' 
i=i 3 1 



N 



T-r sin kj — Ui + iJ 

sin fc,- — ui — iJ 



M 

n 

fe=i 



ui — Uk — 2iJ 
u i — Uk + 2i J' 
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follows that when J > 1/2 the ground state is still the product of rung singlets with 
energy E = —3JL. This is in fact the reference state used in the Bethe ansatz calculation 
and corresponds to the case M = N = for the Bethe ansatz equations. To describe 
an elementary excitation to a spin 1 state we take N=l, M=0 in (H) which yields real 
solutions for the variable k, viz. 



2irr 



0,1,2,...,L-1. 



It is then apparent from the energy expression (§) that for J > 1/2 these elementary 
excitations are gapped. We therefore deduce that J c = 1/2 gives the critical point between 
the gapped and gapless phases of the elementary spin excitations alluded to earlier. 

The model also exhibits elementary bound state excitations which we illustrate in the 
two-site case. For L = 2, N = 2, M = 1 there is a solution of the Bethe ansatz equations 
given by 

u = 0, k\ = —k 2 = arccos(— J) 

which describes an excited state of energy E = — 2J. From the eigenvalue expression 
for (H) we see that this state has zero spin. Such a state has the interpretation of the 
excitation of two bound quasi-particles of opposite spin. 

In order to obtain the solution of this model, we begin with the coupled spin version of 
the Hubbard model as introduced by Shastry [|15j with the imposition of twisted boundary 



conditions. The local Hamiltonian has the form 

h i[i+l) = -<rt%+i) ~ a (i+i) a ^ ~ ^ +r (7+i) - r (t+i) r r 

- | {{at + I){t* + + « m) + J)(rf i+1) + /)) + j 

where {erf, cr?} and {rf, t?} are two commuting sets of Pauli matrices acting on the site 
i. For our convenience an additional applied magnetic field term has been added and an 
overall factor of -1 included. For the twisted boundary term we take 

" j {{al + I)(r[ + /) + K + J) W + I)) + ~- 
The first step is to apply a non-local transformation given by 



i-l 



0{af) = atHr, 



fc=i 



L 



0(rt) = rf II 

k=i+l 



a 
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Under the transformation 9 we yield a new Hamiltonian of the form (^) where the bulk 
two-site operators now read 

hi(i+i) = -°i a (i+i) T i - a (i+i) a i T i ~ ^ r (i+i)°"(i+i) ~~ r (i+i) r i °ft+i) 

- | ((*? + r)(7f + (<rf i+1) + I)(rf m) +/)) + | (6) 

and the boundary term is given by 

L L L L 

hLi = -e-^ai^rt J] rl - e^a+^r! J] r fc ^ - e-^r+rfaj J] o z k - e^r L a z J] 

fc=i fc=i fc=i fc=i 

- \ (K + I)(rl + + K + T) {rl + /)) + j (7) 

An important observation to make is that the boundary term above has non-local 
terms. To accommodate for this, note that we may write 

L 



fc=i 

f[r| = (-1)^ 



fc=i 



where M = J2i-i m ^ N = J2i=i n i and 

m = \{I-° Z ) 

n = I- l -{a z + r z ). 

Since the global operators A4, Af are conserved quantities, we can treat the twisted 
boundary conditions in (|7|) in a sector dependent manner. Letting M and N denote the 
eigenvalues of A4, Af respectively, we now choose 

e i<h = (_ 1 )(iV-M) ) e iH = (_!)M_ (g) 

The validity of making this choice without destroying the solvability stems from the fact 
that states with differing values of M and N are orthogonal independent of the values of 
4>i and 4>2- Hence we may choose different values of </>i and </> 2 for each of the subspaces 
corresponding to a fixed M and N. 

The next step is to now employ a local transformation on the Pauli matrices which 
has the form 

$((7*) = -a + T - - (J-T + + ^((T Z + T Z ) 

$( r z ) = -a+ T - - a-T + - ^(a z + r z ). 



It is worth noting that the above transformation can be expressed 

= TxT' 1 

where 

T = (- + r x -%(-- -L) a z T y - -^=<jv - -^a x r z 

\2 2V2J V 2 2V2J 2V2 2V2 

is a unitary operator. Applying this transformation to the local Hamiltonians (|^,0) gives 
us the local ladder Hamiltonians 

hi( i+ i) = T+a7 +1 + rrat +1 + ((x+r m + a7r+ +1 ) r z a z l+l + j {Bi.fi + a l+1 .f l+1 ) 

and the global Hamiltonian has regular periodic boundary conditions. 

The final step in obtaining (§) is to set J = U/A and perform the transformation 

a, = a t + 1/2 (1 - (-I) 1 ) (n - a,) 
n = Ti + 1/2 (1 - (-1)*) fa - n) 

which has the effect of interchanging the leg spaces on the odd rungs while leaving the 
even numbered rungs unchanged. 

The energy expression for the Shastry model with twisted boundary conditions can 
be obtained through the Bethe ansatz. The result is 

N 



such that the kj satisfy the Bethe ansatz equations 

(-l)^e-^exp(zLA;,) = - f[ Sin k J - * + iJJ jf J = 1,2,...,A 

sin kj — m — lu /A 

yr sinkj-m + iU/A = 1 sN c -i(*x-* a )T\ u i- u k-iU/2 , = 12 M 
J-l sin kj - ui - iU/A v ; Ai uj - u k + iU/2 J ' 

j — 1 ■ k — 1 

An important point here is that the numbers M and iV above have precisely the 
same meaning as the interpretation presented earlier; i.e. they are the eigenvalues of the 
conserved operators M. and N '. Consequently, we need only substitute the values of eq. 
(§) into the above energy expression and Bethe ansatz equations which gives us with 
the parametrization J = U/A. 
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